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Abstract 

The S-matrix of diffractive scattering is diagonalized in terms of the colour dipole-dipole cross 
section. Recently, we have shown that the dipole cross section satisfies the generalized BFKL 
equation. In this paper we discuss the spectrum and solutions of our generalized BFKL equation 
with allowance for the finite gluon correlation radius R c . The latter is introduced in the gauge 
invariant manner. We present estimates of the intercept of the pomeron and find the asymptotic 
form of the dipole cross section. We consider the difference between the BFKL and GLDAP 
evolutions and conclude that the GLDAP evolution remains a viable description of deep inelastic 
scattering well beyond the kinematical range of the HERA experiments. We suggest methods of 
measuring the pomeron intercept in the HERA experiments. 
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1 Introduction. 



The asymptotic behavior of high-energy scattering in perturbative QCD is usually discussed in 
terms of the Balitskii-Fadin-Kuraev-Lipatov (BFKL) pomeron [1-3]. The BFKL equation is the 
integral equation for the differential density of gluons, so that formulation of predictions for the 
observable quantities like total cross section is not straightforward. Furthermore, the original 
BFKL equation is the scaling equation, and when discussing the realistic QCD one encounters 
the difficult task of introducing the infrared cutoff in the BFKL equation. The scaling BFKL 
equation has a continuous spectrum and corresponds to the vacuum singularity described by a 
cut in the complex angular momentum plane. One of the pressing issues in the theory of the 
perturbative QCD pomeron is whether the infrared cutoff and the running QCD coupling can 
change the spectrum of the pomeron and produce isolated poles. 

In this paper we discuss the spectrum and solutions of the generalized BFKL equation directly 
for the total cross section, derived by us recently [4,5]. Our equation allows to introduce a finite 
correlation radius for the perturbative gluons in the gauge invariant manner. The starting point 
of its derivation is the technique of multiparton lightcone wave functions, developed by two of the 
present authors [4,6]. The principal observation is that the transverse separations /% = f>i — Pj 
and the lightcone momentum partitions Zi of partons in the many-body Fock state are conserved 
in the scattering process. Interaction of the (n+2)-parton Fock state is described by the lightcone 
wave function ^ n+2 (p n +2, z n+2 , ...,p\,Zi) and the (n + 2)-parton cross section cr n+2 (p n +2, ■ ■■,Pi), 
which are perturbatively calculable [4] (here n referes to the number of gluons in the Fock state). 

To the lowest order in the perturbative QCD one starts with the qq Fock states of mesons 
(qqq state for the baryons) and with the scattering of the (A) projectile and (B) target colour 
dipoles of transverse size ta and r B (here Fa, r B are the two-dimensional vectors in the impact 
parameter plane) 



OQ(r A ,r B ) = — 




a 2 s 1 — exp(— ikrA) 1 — exp(ifcre) 



(1) 



Here 





(2) 
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is the running strong coupling, /3q = 11— ~Nf = 9 for Nf = 3 active flavours, and in the integrand 
of the dipole-dipole cross section oil must be understood as as(max{fc 2 , ^})as(max{/c 2 , ^r}), 
where C ~ 1.5 [6]. Introduction of the infrared freezing of the running coupling and of the 
correlation radius for gluons R c — 1/^g is discussed below. In terms of the dipole-dipole cross 
section ([I]) the perturbative part of the total cross section for the interaction of mesons A and 
B equals 

a {pt \AB) = ({cr(f A ,r B )) A ) B = J dz A d 2 f A dz B d 2 f B \^(z A ,f A )\ 2 \^(z B ,f B )\ 2 a (f A ,f B ) . (3) 

The irrefutable advantage of the representation (|3p is that it makes full use of the exact di- 
agonalization of the scattering matrix in the dipole-size representation. (Hereafter we discuss 
ao(r,R) averaged over the relative orientation of dipoles, as it appears in Eq. (|3|).) Notice the 
beam-target symmetry of the dipole-dipole cross section, 

cro(r A , r B ) = <J (r B , r A ) , (4) 

and of the representation (|3]). 

The increase of the perturbative component of the total cross section comes from the rising 
muliplicity of perturbative gluons in hadrons, n g oc logs, times Aa g - the change in the dipole 
cross section for the presence of gluons: Aa^ ~ n 9 A<r 9 . (Here s is a square of the cm. s. energy.) 
The n-gluon Fock components of the meson give contributions oc log n s to the total cross section. 
The crucial observation is that the effects of higher order Fock states of the interacting hadrons 
can be reabsorbed into the energy dependent dipole-dipole cross section with retention of the 
representation @. This energy-dependent dipole cross section satisfies the generalized BFKL 
equation derived in [4,5], investigation of properties of which is the subject of the present paper. 
The presentation is organized as follows: 

In section 2 we briefly review the derivation of our generalized BFKL equation for the dipole 
cross and in section 3 discuss its BFKL scaling limit. In section 3 we also comment on the 
introduction of the infrared cutoff into the scaling BFKL equation. The subject of section 
4 is the conventional QCD evolution as the limiting case of our generalized BFKL equation. 
The impact of the finite correlation length for gluons and of the running QCD coupling on the 
spectrum of the generalized BFKL equation is discussed in detail in section 5. Our principal 
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conclusion is that the generalized BKL kernel has a continuous spectrum, so that the partial 
waves of the scattering amplitude have a cut in the complex angular momentum plane. We also 
comment on the conditions under which the pomeron can have a discret spectrum. In section 6 
we find the form of the dipole cross section for the rightmost singularity in the j-plane and its 
intercept as a function of the gluon correlation radius R c . The subject of section 7 is a transition 
from the GLDAP evolution to the BFKL evolution with increasing energy. The beam-target 
symmetry of the dipole-dipole cross section and the admissible form of the boundary conditions 
for the generalized BFKL equation are discussed in section 8. In section 9 we comment on the 
restoration of the factorization property of the asymptotic cross section. In section 10 we suggest 
practical methods of measuring the intercept of the pomeron in the HERA experiments on deep 
inelastic scattering. In the Conclusions section we summarize our basic results. 



2 Generalized BFKL equation for total cross section 

Now we sketch the derivation of our generalized BFKL equation [4,5] for total cross sections. 
Unless specified otherwise, we consider the contribution to the hadronic scattering from the ex- 
change by the perturbative gluons and suppress the superscript (pt). The perturbative qqg Fock 
state generated radiatively from the parent colour-singlet qq state of size r has the interaction 
cross section [4] 

9 1 

^(r, Pl ,p 2 ) = -[<7 (pi) + cr (p 2 )} ~ g<7 (r) , (5) 

where p ij2 are separations of the gluon from the quark and antiquark respectively, p 2 = P\ — r. 
(Hereafter we suppress the target variable Tb and for the sake of brevity use r = ta-) The 
cross section cr 3 (r, p±, p 2 ) has the gauge invariance properties of (Jz{r, 0, r) = <Jz{r, r, 0) = o"o(r) 
and cr 3 (0,p, p) = §o" (p) [4]. The former shows that when the gluon is sitting on top of the 
(anti) quark, the qg(qg) system is indistinguishable from the (anti) quark. The latter shows that 
the colour-octet qq system of the vanishing size is indistinguishable from the gluon, and | is 
the familiar ratio of the octet and triplet couplings. The increase of the cross section for the 
presence of gluons equals 

9 

Aa g (r, Pu p 2 ) = a 3 (r, Pl , p 2 ) - a (r) = -[a ( Pl ) + a (p 2 ) - a (r)} , (6) 
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The lightcone density of soft, z 9 < 1, gluons in the qqg state derived in [4] equals 



\$i{r, Pi, p 2 , z g )\ 2 = —t^Pg 



^ / (min)\ I \ Pi / (min)\ Ty - / \P2 

g s [r\ ^K^jiGpx) g s (r 2 ')K 1 (p G p 2 ) — 



(7) 

Pi P2 

Here <7s(r) is the running colour charge, as(r) = gf(r)/47r is the running strong coupling, 
^(rrmi) _ mm { rj ^ ^ Ki(x) is the modified Bessel function, z g is a fraction of the (lightcone) 
momentum of the qq pair carried by the gluon, and / dz g /z g = log(s/so) = £ is the parameter 
of the Leading-log^) approximation and gives the usual logarithmic multiplicity of radiative 
gluons. The wave function (|7p counts only the physical, transverse, gluons [4]. 
If n g (r) is the number of perturbative gluons in the dipole r, 



n g {r) 



J dz g d 2 pi |$i(f, pi,p 2 ,z g )\ 2 , (8) 



then the weight of the radiationless qq component will be renormalized by the factor [1— n g (z, r)]. 
As a result, the total cross section with allowance for the perturbative gluons in the beam dipole 
A takes the form 



<7(€>r) = [l-n g (z,r)]ao(r) + J dz g d 2 p |$i(f, p u p 2 , z g )\ 2 a 3 (r, p u p 2 ) 
= (T (r) + J dz g d 2 p! \<$> 1 {f 1 p u p 2} z g )\ 2 /\ag{r,p u p 2 ) =<7 (r) + <7i(r)f = [1 + £/C®]<7 (r) , (9) 

where the kernel K is defined by [4,5] 



01 (r) = K, ® a (r) 



3 



in 3 



J d 2 pi p 2 G 



, (min) \js I \ Pi / (min) \js-l \ P% 

g s (r\ ; )^i(PgPi) 9s(r 2 'jK^pcp^ — 



Mpi) + cr {p 2 ) - a {r)} . (10) 



Pi P2 

Eq. (|9p shows that the effect of gluons can be reabsorbed into the generalized, energy-dependent, 
dipole cross section cr(£,r). To higher orders in £, this generalized dipole cross section can be 
expanded as <r(£, r) = J2 n =o ^[ cr «-( r )£ , \ where <J n+ i — K ® cr n , so that 

^■ = K9^,r) (11) 

is our generalization of the BFKL equation for the dipole cross section. 

The colour gauge invariance of the presented formalism is noteworthy. Firstly, the dipole- 
dipole cross section coO^a, re) vanishes at ta — > or r# — > 0, because gluons decouple from 
the colour-singlet state of vanishing size. Secondly, for the same reason the wave function (|7|) 



vanishes at r — > 0. Thirdly, Aa g (r, pi, p?) — > when p\ —>■ (or p 2 — > 0), since by colour 
charge conservation the quark-gluon system with the gluon sitting on top of the (anti) quark is 
indistinguishable from the (anti) quark, and the interaction properties of such a qqg state are 
identical to that of the qq state. Therefore, our introduction of the finite correlation radius 
for gluons R c = l/pc, which takes care of the perturbatie gluons not propagating beyond the 
correlation radius R c , is perfectly consistent with the gauge invariance. For instance, Eq. ([T0|) 
supports the gauge invariance constraint a n (r) —>■ at r — > 0, to all orders n. (We do not have 
a complete proof of the gauge invariance, though. Neither can we prove that our prescription is 
unique.) 

The renormalization of the weight of the radiationless qq Fock state in Eq. ([]) in a simple and 
intuitively appealing form takes care of the virtual radiative corrections (in the BFKL formalism 
[1,2] these very radiative corrections are responsible for the reggeization of gluons). Notice, that 
although the multiplicity of radiative gluons n g (r), which appeared in the intermediate stage of 
the derivation of our Eq. (|T0l) , is formally divergent, the generalized dipole cross section a(£,r) 



and Eqs. pOL 1 1 ) are both ultraviolet and infrared finite. 



3 The BFKL scaling limit 



In the BFKL scaling limit of r, pi, p 2 "C R c and of the fixed as 



/4 



K^pGpi)— - K 1 (p G p 2 ) — 

Pi P2 



1 1 ' 

pipi 



;i2) 



the kernel K, becomes independent of the gluon correlation radius R c and with the fixed as it 
takes on the scale-invariant form. The corresponding eigenf unctions of Eq. (|II]) are 

E{uo,i,r) = (r 2 )^ +<JJ exp[£A(a/)] 



(13) 



with the eigenvalue (intercept) [here r = rn, p\ = rx and p 2 = r(x + n)} 



AM 



lim^f / d 2 x 

e^O 27T 2 



2{x 2 ) 



[x 2 + e 2 ] [(x + n) 2 + e 2 ] n Jo 



3a s f 1 , Z2~ 
az — 



+ Z2 +LJ - 2z 



z(l - z) 



3a s 
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[2*(1) 



-u) 



(14) 



Here \l/(x) is the digamma function, and we have indicated the regularization which preserves 
the symmetry of the kernel /C. The final result for A(o;) coincides with eigenvalues of the BFKL 
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equation found in [1-3]. The solution ( pT3|) corresponds to the singularity in the complex- j plane, 
located at j = 1 + A (a;). The rightmost singularity is located at j = ajp = 1 + Ajp, where 

AjP = A(0) = m^a s . (15) 

71 

When uj is real and varies from — | to and to |, also the intercept A (a;) is real and varies 
from +oo down to A(0) = Ajp and back to +oo, along the cut from j = 1 + Ajp to +oo in the 
complex angular momentum j plane. If ou — iv and v varies from — oo to and to +oo, then 
the intercept A(iz/) is again real and varies from — oo up to A(0) = An> and back to — oo, along 
the cut from j = — oo to j = 1 + Ajp in the complex j-plane. The choice of the latter cut is 
appropriate for the Regge asymptotics at £ ^> 1 and the counterpart of the conventional Mellin 
representation is 

/+oo r+oo 
dv f{v)E{iv,r,£) = r I duf(u) exp[2wlog(r)] exp(A(w)£) . (16) 
-oo J — oo 



The Mellin transform is obtained from Eq. ( |16D if one changes the integration variable from the 
v to the angular momentum j = 1 + A(iu). The spectral amplitude f(u) is determined by the 
boundary condition cr(£ = 0,r): 

f{y) = I f dr ^fl exv [-2ip\og(r)] . (17) 

Eqs. ( P^|JT7| ) give a nontrivial connection between the r-dependence of the total cross section 
and its energy dependence. In the BFKL scaling regime, the rightmost j-plane singularity 
corresponds to the asymptotic dipole cross section crip(£,r) = <xip(r) exp(£Aip), where 

crip(r) oc r . (18) 

More direct correspondance between our equations (]10|,[TTD in the scaling limit of fie — > and 



the conventional BFKL equation can be established if one rewrites Eqs. (|10|JTl| ) as an equation 
for the function 

9 (M = ^4 (19) 

which in the BFKL scaling limit is simply the density of gluons g(x, k 2 ) at the Bjorken variable 
x = xoexp(— £) and the virtuality k 2 ~ 1/r 2 , where xq ~0. 1-0.01 corresponds to the onset of the 
leading-log(^) approximation (the relation ( |I9"D holds also for the running strong coupling as = 
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J d 2 k f_ [g(£, k 2 ) + gfc (p- kf) - g&p 2 )} . (20) 



as(r), for the detailed derivation see [4,7]). The transformation from the size-r* representation to 

the momentum-p representation can easily be performed making use of the conformal symmetry 

of the BFKL equation [3] . Namely, by virtue of the conformal symmetry r can easily be traded 

for p and cr(£,r) can be traded for g(£,,p 2 ) with the result 

dg(£,p 2 ) _ 3a s r p_ 2 

d£ " 2tt 2 J k 2 (p- k) 2 

The factor p 2 in the integrand of fl2Tj) is crucial for supporting the gauge invariance constraint 
9(£,,P 2 ) — > at p 2 ^ 0, which is the counterpart of the gauge invariance constraint <x(£, r) — > 
at r — > in the ^-representation. If one introduces </>(£, k 2 ) = g(£, k 2 )/k 2 , then the original form 
of the BFKL equation will be recovered from Eq. (|20"D (in the scaling limit 0(£, k 2 ) satisfies the 
same equation as the more often considered dg(£, k 2 )/dk 2 ). 

Introduction of the infrared regularization into the BFKL scaling equation (|20| ) is one of the 
pressing issues in the theory of the perturbative QCD pomeron. We are going to use neither 
Eq. pD] ) nor any of its modifications to be discussed below; we resort to our Eqs. (|TI],[n]), which 
allow to treat effects of both the gluon correlation radius R c and the running strong coupling 
in the consistent and manifestly gauge invariant manner. Nonetheless, in order to establish a 
connection with the existing literature, we comment on the infrared regularization of the BFKL 
equation fl2~0|) in the momentum representation. 

Our form (p0|) of the BFKL equation is more convenient for these purposes, as it allows the 
clearcut separation of propagators 1/k 2 and l/(p — k) 2 from the factor p 2 in the nominator in 
Eq. (20). The latter, as well as the factor r 2 in Eq. ([12]) , has its origin in cancellations of the 



radiation of soft gluons by the colour singlet states. Therefore, the infrared regularization of 
Eq. (p0|) which plausibly respects the gauge invariance, is the substitution 

^ - V 2 (21) 



k 2 (p — k) 2 [k 2 + in 2 ] \{p — k) 2 + m 2 j 
which treats the both gluon propagators on the equal footing. In the scaling limit of p 2 ^> m 2 , 

the function ipi^P 2 ) = 9{£,iP 2 )/{p 2 + m2 ) is similar to, and satisfies the same equation as, 
dg(£,p 2 )/dp 2 . The equation for ip(£,,p 2 ), which follows from (|20| ) subject to the substitution 
(|21|), reads 

^,k 2 ) + m(p-k) 2 ) (p 2 + m 2 )^,p 2 



dipj^p 2 ) = 3o5 _ P 2 i (flj{ 
<9£ 27r 2 p 2 + m 2 



(p — k) 2 + m 2 k 2 + m 2 [(p — k) 2 + m 2 ] [k 2 + m 2 



71 



2 p 2 + m 2 J L(p-A;) 2 + m 2 [(p - k) 2 + m 2 } [k 2 + (p — k) 2 + 2m 2 ] 



(22) 



The asymmetric form in the second line of Eq. ( p2|) is obtained from the symmetric form in the 
first line of Eq. (|22| ) if one uses the identity [2] 

[ ^ -of (23 ) 

J [k 2 + m 2 ][(p- k) 2 + m 2 ] J [k 2 + m 2 }[(p-k) 2 + k 2 + 2m 2 } 

(Erroneous form of this identity is cited in [8]). The asymmetric form of Eq. (|22|) at m 2 = 
is precisely the Balitskii-Lipatov form [2] of the BFKL equation, which is the one customarily 
used in the recent literature. 

One often imposes on this asymmetric version of the BFKL equation the infrared cutoff 
k 2 > k 2 [9-11], which is illegitimate. Indeed, such a cutoff does not ensure, and as matter of fact 
manifestly breaks, the condition of the identical infrared cutoff of different gluon propagators, 
as it breaks the initial symmetry property of the BFKL equation ( p0|) (for a similar criticism 
of the Collins-Kwiecinski (CK) sharp cutoff k 2 > k 2 [9] see also Collins and Landshoff [10]). 
Furthermore, such a cutoff is evidently suspect from the point of view of the gauge invariance. 
One must rather use the equation ([22|) which introduces the infrared regularization with retention 
of the symmetric cutoff of the gluon propagators and in the manifestly gauge invariant manner. 
Balitskii and Lipatov [2] did not encounter these problems, as they introduce the massive vector 
mesons in a consistent manner, using the spontaneous breaking of gauge symmetry. Introduction 
of massive gluons into the BFKL equation was also discussed by Ross and Hancock [11]. The 
Ross-Hancock prescription is different from ours in Eq. (^), manifestly breaks the symmetry of 
the cutoff of gluon propagators and, in contrast to our Eq. (P2"|), does not support the p 2 — ► 
gauge invariance constraints. 



4 The GLDAP limit 

The limit of large but finite £ and very small r 2 R 2 is often referred to as the Double- 
Leading-Logarithmic Approximation (DLL A). Here R is either the target size, R ~ re or the 
correlation radius R c if R c < tb- This is the limit typical of the deep inelastic scattering, 
and it is usually described by the GLDAP evolution equations [12]. Introduction of the gluon 
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correlation radius R c and of the running QCD coupling poses no problems in the DLLA. In the 
DLLA, £ = log(— ), where x is the Bjorken variable, and xq < 0.1 — 0.01 is the threshold for 
the Leading-log (^) approximation. The BFKL equation contains the GLDAP equation as the 
limiting case. In the DLLA the kernel /C takes a particularly simple form [4] 



t n ^ / \ 3r 2 a s (r) f R2 d 2 p . . . A . 
a n+1 (r) =K® a n (r) = ^ / -f cr n (p) , (24) 

7T J r 2 p 



which is equivalent to the GLDAP evolution equation. The dominant contribution comes from 
the DLLA ordering of sizes r 2 <C p 2 <C R 2 , which justifies factoring out «s(r) in Eq. (|24l) . In 
the same DLLA the boundary condition, the dipole-dipole cross section (|3]), equals 

<r(f = 0, r) = <7 (r, i?) « Cr 2 a 5 (r)L( J R, r) , (25) 

where L(R,r) is the parameter of the Leading- log (r 2 ) approximation (the Leading-log(Q 2 ) ap- 
proximation in the more conventional momentum representation), 

~a s (RY 



L(R,r) « log 

[a s {r) 

Iterations of Eq. fl24|) subject to the boundary condition fl25|) give the DLLA solution [4] 



(26) 



a n+1 (r) = /C ® tr n (r) = — — ■ -^^(i?, r)a„(r) . (27) 

71 + 1 A) 

The corresponding generalized dipole cross section has the DLLA asymptotics 

°dlla{Z, r) = a(0, r) £ ~ r 2 a 5 (r)L(r, , ( 28 ) 

where 7/ is the expansion parameter of DLLA, which is a product of the leading-log (s) and the 
leading-log(r 2 ) expansion parameters: 

V = ^L(r,R) = UL(r,R). (29) 

Notice, that the DLLA solution implicitly invovles an assumption that at large r ~ R, the cross 
section (Tdlla(£, r ) is approximately constant. 

The large-n behavior of the DLLA iterations (]27 ) must be compared to 



cr n +i(r) = A-spa n (r) (30) 
10 



for the leading BFKL solution, which suggests the criterion for the breaking of the DLL A and/or 
GLDAP evolution: At large r\ the DLLA cross section ( p8|) is dominated by the contribution 
from large n ~ The DLLA to GLDAP evolution breaks when 



The pattern of this DLLA breaking in the case of the running as(r) will be discussed in more 
detail below. Here we only wish to emphasize that in the physically interesting case of the 
running coupling and of finite correlation radius for the perturbative gluons R c the intercept 
Aip is the constant, which does not depend on r. Therefore, the boundary line x = x c (r) in the 
(r, £ ) plane between the BFKL evolution and the GLDAP evolution to DLLA is given by 



It is different from the much discussed erroneous boundary suggested in [8]. We shall comment 
more on the transition from the DLLA to GLDAP evolution to the BFKL evolution below, now 
we concentrate on the investigation of the spectrum and eigenfunctions of our generalized BFKL 
equation. 

5 The spectrum and eigenfunctions of the generalized 
BFKL equation 

5.1 The Green's function, diffusion and the 'Schrodinger' equation 

The generic solution of Eq. (|TTD can be written as 



L(R,r) _ L(R,r) < 3 
n y/rj ~ 4 



ip • 



(31) 




(32) 




(33) 



where in the BFKL regime the evolution kernel (Green's function) equals 




(34) 



Since A(iv) has the maximum at v — 0, 



A(zi/) = 



(35) 
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the large-£ behavior of the Green's function can be evaluated in the saddle-point approximation: 



v(c /\ exp(Aip^) 
K (£ , r, r ) oc -j= exp 



(logr 2 — log(r^ 2 ,2 



(36) 



2£A"(0) 

Since the BFKL eigenfunctions (|T3| ) are oscillating functions of r, a priori it is not obvious that 
an arbitrary solution ( |I6"D will be positive-valued at all £ and r. The BFKL kernel K(^,r,r') of 
Eq. (|34]) is manifestly positive- valued at large £, so that starting with the positive- valued <t(0, r) 
one obtaines the positive- valued asymptotical cross section cr(£,r). 

The 'diffusion' kernel (|3~4D makes it obvious that the BFKL scaling approximation is not 
self-consistent in the realm of realistic QCD: starting with <r(£ = 0,r) which was concentrated 
at the perturbative small r <R -C R c one ends up at large £ with cr(£,r) which extends up to 



the nonperturbative r ~ i?exp(y'^A"(0)) > R c . Therefore, introduction of a certain infrared 
regularization in the form of a finite gluon correlation length R c is inevitable, and hereafter 
we concentrate on effects of finite R c on the intercept Ajp. Interpretation of fie — 1/Rc as 
an effective mass of the gluon suggests the infrared freezing of the strong coupling as{r) at 
r > Rf ~ R c . In the numerical analysis we use the running coupling 

Att 

a s( r ) = 7 " — \~ ' (37) 



/5olog 



c 2 

^QCD r 



where C ~ 1.5 [6]. At large r we impose the simplest freezing cts(r > Rf) = ctg = 0.8. This 
corresponds to the freezing radius Rf ~ 0.42f and/or the freezing momentum kf w 0.72GeV/c 
in the momentum representation (fj). The more sophisticated smooth freezing can easily be 
considered, but it will be obvious that our principal conclusions do not depend on the form 
of the freezing. The gluon correlation radius R c and the freezing coupling provide the minimal 
infrared regularization of the perturbation theory, and in the sequel we consider the so-regularized 
generalized BFKL equation ( |10] , []~T1) as applicable at both small and large radii r. 

Although only the case of R c ~ Rf is of the physical interest, a study of (albeit unphysical) 
limiting cases Rf — > at finite R c , and of finite Rf at R c —>■ oo, is instructive for the insight into 
how the spectrum of the j-plane singularities is modified by the infrared regularization. The 
corresponding analysis is greatly facilitated by the observation, that in view of Eqs. Q3"5|,|3T)D the 



large-£ behavior of solutions of the BFKL limit of Eq. flTT|) is similar to the large-£ behavior of 
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solutions of the 'Schrodinger' equation 

A"(0) d 2 



+ V(z) 



5$ 



(38) 



2 <9z 2 

for a particle of mass M = 1/A"(0) in the potential V(z) = — A(0). Here the coordinate 
z = logr 2 , the 'wavefunction' <&(z,£) = cr(£,r)/r and the intercept equals the 'energy' e taken 
with the minus sign. 



5.2 Finite correlation radius R c , fixed coupling as 

The first limiting case of Rf —>■ corresponds to introduction of a finite gluon correlation radius 
R c at fixed as- In this limiting case the intercept Ajp will still be given by the BFKL formula 
fli~5|) . Indeed, on the infinite semiaxis logr < logi? c the kernel K, retains its scaling properties, 
the corresponding eigenfunctions will be essentially identical to the set (|T3"D, the spectrum of 
eigenvalues will be continuous and the cut in the j-plane will be the same as at R c —>■ oo. This is 
particularly obvious from the Schrodinger Eq. (p8|), since on the the semiaxis z < z c = \ogR 2 c the 
potential V(z) = — Ajp is flat, and the corresponding Schrodinger operator (|38|) has a continuous 
spectrum starting with the minimal energy e = — Ajp. Evidently, this property of the spectrum 
of the Schrodinger operator (|38| ) does not depend on the details of how the gluon correlation 
length R c is introduced. This is a reason why the similar conclusion on the spectrum of the 
infrared-cutoff BFKL equation was reached in [9] in a model with the above criticized sharp 
k 2 > k 2 cutoff in the momentum representation. 

The behaviour of solutions at large r 3> R c requires special investigation. In this region 
the term cx K 1 (fioPi)K 1 (fi G p 2 ) in the kernel K is exponentially small, which is related to the 
exponential decay of the correlation function (the propagator) of perturbative gluons. (Such an 
exponential decay of the gluon correlation function with the correlation radius R c = 0.2 — 0.4f 
is suggested by the lattice QCD studies, for the review see [13]. We remind that Ki(x) in the 
kernel K, comes from the gradient of the gluon correlation function oc K (x)). Then, at large r 
the kernel /C of our Eq. (0) will be dominated by the contributions from pi < R c p 2 ~ r and 
from p 2 ^ R c "C pi ~ r, will have the limiting value which does not depend on r, and Eq. (|TT1) 
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takes the form 

o n +i{r) = -^f J d 2 pi /4^i(^gPi) 2 K(Pi) + ^nifto) - °V>( r )] • ( 39 ) 
This equation has a continuum of solutions with the large-r behavior of the form 

E + ((3, £, r) = [a{(3) + exp(z/^ G r)] exp[5(/?)f ] . (40) 

(Apparently, one can put a(/3) = 0, but this is unimportant for the subsequent analysis). 

When continued to small r through the region of r ~ R c , the plane waves in the linear-r 
space transform into the plane waves in the logr 2 -space (times the overall factor r), E(ii>, £, r) = 
r exp(iz/log r 2 ) of Eq. flT3|), so that S(j3) = A(iv). Evidently, the two real and the node-free 
solutions with u = and (3 = must match each other, so that the rightmost j-plane singularity 
with the intercept A#> Eq. ( |I~5|) must correspond to 0ip(r) oc r at r <C R c and o-jp(r) oc const at 
r ^> R c . We wish to emphasize that although such a 0ip(r) extends to large r, Eq. (|39D makes it 
obvious that the intercept Ajp is controlled by the behavior of ajp(r) at r ~ i? c and is relatively 
insensitive to the large-r behavior of ajp(r). 

The intercept 8((3) will have a maximum at /3 = 0. Then, repeating the derivation of the 
Green's function ( JT7\ ) one can easily show that at large r, r' the Green's function of Eq. (^) has 
the behavior 

exp(AjpO / u 2 c (r-r') 2 \ 

A («-"' K ^V^ ex H"^(or)' (41) 

which is reminiscent of the familiar multiperipheral diffusion. Indeed, at large r ^> 1/ u G & sort 
of the additive quark model is recovered, in which the (anti) quark of the dipole developes its 
own perturbative gluonic cloud, and the quark-quark scattering will be described by the multi- 
peripheral exchange of massive vector mesons. The Green's function (|28|) shows the emergence 
of the Regge growth of the interaction radius in this limit. It also shows that the large-£ behavior 
of solutions of our generalzed BFKL equation in this region will be similar to solutions of the 
Schrodinger equation 



5"{0) d 2 



U{r) 



d 

$ = -— $ = e$ (42) 



2fiQ dr 2 

with the potential U(r) = —6(0) which is flat at r > R c . 

The scaling limit of Eq. (|TT|) does not have the localized solutions, see the eigenf unctions (p!3|). 
An important observation is that also Eq. (^) does not have the node-free localized solutions, 
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which vanish at large r. Indeed, if such a solution had existed, then the l.h.s. of Eq. ( |39|) would 
have vanished at r — > oo. For the same reason, we can neglect a(r) and cr{p2) in the integrand of 
the r.h.s. of (^) and will be left with the positive-valued, and r-independent, contribution from 
er(pi). This shows, that (f40| ) gives a complete set of solutions. Consequently, the so generalized 
BFKL equation only has the continuous spectrum, and the partial-wave amplitudes only have 
the branching point singularity (cut) in the complex- j plane. 

We wish to emphasize that none of the above conclusions changes, if Ki(fiar) is changed 
for any function which vanishes sufficiently steeply at large r, the only condition is that it must 
have the short-distance 1/r behavior dictated by the 1/k 2 ultraviolet behavior of the perturbative 
gluon propagator. 



5.3 The running and freezing coupling, massless gluons 

The second interesting case is of R c — > oo, i.e., the case of massless gluons, \iq — > 0, but with 
the running coupling which freezes, a§(r) = org = atg(Rf), at finite r > Rf. In this case the 
scaling invariance of the kernel K. is restored on the infinite semiaxis logr > logi? c , where the 
eigenf unctions E+(iu, £, r) are again essentially identical to the BFKL set fli~3f) , the spectrum of 
eigenvalues will evidently be continuous, the intercept A#> will be given by the BFKL formula 
QT5D with as = org and the partial waves of the scattering amplitude will have the cut in the j- 
plane identical to that for the original BFKL equation. Indeed, the corresponding 'Schrodinger' 
equation has the continuous spectrum starting from e = — A(0), and this conclusion is 
evidently insensitive to the exact shape of the potential V{z) at z < Zf. Here we agree with 
Ross et al. [11,14] and disagree with Lipatov [3], who concluded that introduction of the running 
coupling is by itself sufficient for transformation of the continuous spectrum of the pomeron to 
the discret one. Apparently, the origin of this conclusion is in that Lipatov [3] restricted himself 
to a very special form of the solution at large r, rather than using the continuum solutions (flOl). 

The continuation of the BFKL solutions ( |i~3"D to the region of z < Zf = logR 2 poses no 
problems. Let us start with the quasiclassical situation when < 1, so that ots{z) = 

cus(Rf) [1 + (Po/4:Tr)a!g(Rf)(zf — z)}~ 1 is a slow function of z. In this case the slowly varying 
running coupling can be factored out from the integrand of Eqs. ( |10yTT| ). This suggests that to 
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a crude approximation one may make a substitution of the fixed coupling as in Eq. (0) for the 
running coupling cts(r). Consequently, Eq.(^) takes the form of the 'Schrodinger' equation for 
a particle with the slowly varying mass 

~ ^ ■ (43) 
1 ; A"(0) a s (r) K } 

in the slowly varying and monotonically decreasing in the magnitude potential 

^~- a *;t7^' (44) 

a s {Kf) 

which is flat, V(z) = — Ajp, at z > Zf. Evidently, the solutions E+(iu,£,r) with A(iu) > 
, i.e., with the negative energy e, will have the underbarrier decrease at z — > — oo, whereas 
the solutions with A(iv) < will be continued to z < as plane waves. This suggests that 
the eigenfunction <Jjp(r) for the rightmost singularity will decrease at r — > faster than the 
oc r 1 solution ([TBI) for the fixed-as scaling BFKL regime. The case of large values of the frozen 
coupling ag must be qualitatively the same. 

In section 4 we have discussed the conventional DLLA solution of the GLDAP equation, 
which has a limited applicability domain. Besides this conventional DLLA solution, in the weak 
coupling limit of as(r) 1 the GLDAP equation (|24]) has the new one-parametric family of 
solutions which have the Regge behavior oc exp(A£). It can best be described in terms of the 
density of gluons 



xg(x, r) 



MM 



(45) 



7r 2 r 2 cts(r) 

After factoring out as{ r ) from the kernel K, which we shall justify a posteriori, our generalized 



BFKL equation (TTTJ) takes the form 



dG{i,r) 



d 2 px 



as(Pi)G(C,Pi) , 0! S (p 2 )G(£,p 2 ) r 2 a s {r)G(£,r) 



Pi 



Pi 



P\p\ 



(46) 



A posteriori, we can show that the last term oc as(r)G(^,r) can be neglected at a$(r) 1 
Then, Eq. ([|7]) takes the form 



4 rHR,r) 
3 Jo 



dL{R,p) G{ti,p) 



which has the one-parametric family of the small-r, large-L(i?, r) eigenfunctions 



e-(7>^>0 = exp[7L(i?,r)]exp(A0 



a sir 



exp(AO 



(47) 



(48) 
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where the exponent 7 is related to the intercept A by 



7=4- (49) 



The corresponding dipole cross section equals 



cr_(7, r) = r 



7-1 

2 



exp(A£) , (50) 



_a s (r)_ 

In contrast to the conventional DLLA solution ( P%|) which summs the leading powers [L(R, r)£] n , 
our new solutions P8|j50| ) do manifestly summ all powers L(R, r) k . The accuracy of this solution 



can be understood making use of (f48|) in Eq. (pE6|): the neglected term oc as(r)G(£,r) gives the 
oc as(r) correction to the solution ([48|). 

The small-r considerations alone do not fix the intercept A and the exponent 7, they are 
determined from the matching the solution e_(7, £, r) with the large-r solution E + (ii>,!;,r). 

5.4 The strong coupling freezes at Rj ~ R c . 

This realistic case of the major interest is a combination of the two previous cases. It is convenient 
to start from the region of r > R c , Rf. The only change from Eq. (^) will be that the running 
coupling as must be absorbed into the integrand 

3 r 

a n+1 (r) = — J d 2 pi as{pi)p 2 G Ki(p G pi) 2 [a n (pi) + a n (p 2 ) - cr n (r)] , (51) 

This equation has a continuum of solutions of the form ( [40]) with the spectral parameter — 00 < 
(3 < +00, and its large-r Green's function and the corresponding Schrodinger operator will be 
similar to ( |4TD and (|42D, respectively. Therefore, Eq. (|51|) has the continuous spectrum and 
the scattering amplitude will have a cut in the complex j-plane. Continuation of the large- 
r solutions to small r is not any different from that discussed in section 5.3. Evidently, the 
intercept Ajp of the rightmost singularity in the j-plane corresponds to the eigenvalue 8{j3) with 
(3 = 0, Ajp = 6(0), as the corresponding eigenfunction does not have a node. For this solution 
&Tp{£,,r) = const vs. r at large r, and decreases faster than oc r 1 at r — » 0. From Eq. fl5T|) it 
is obvious, that the intercept Aip is mostly controlled by the contribution from r ~ R c . Also, 
repeating the considerations of section 5.1, we can exclude the node-free localized solutions of 



Eq. ( |5TD with er(r) vanishing at large r. This shows, that our generalized BFKL equation only 
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has the continuous spectrum and the partial-wave amplitudes only have the cut in the complex 
j-plane. 



5.5 Is the discret spectrum of the pomeron possible? 

The discret spectrum of the pomeron, i.e., generation of isolated poles in the complex- j plane, 
comes along with the existence of localized solutions of the generalized BFKL equation. The 
above elimination of localized solutions is rigorous in the framework of our minimal infrared 
regularization. Technically, it is based on the kernel K being finite at r — > oo, see Eqs. ([39],[5l]), 
which gives rise to a(r) ~ const at r > R c and to the continuum of solutions of the form (flOl). 
The interaction picture which emerges at r > R c has much intuitive appeal: each well separated 
quark of the beam (target) dipole developes the perturbative gluonic cloud of its own and a sort 
of the additive quark model is recovered. At first sight, the localized solutions with the vanishing 
perturbative total cross section for colour dipoles of large size, cr(r) — > at r — > oo, look quite 
unphysical. 

Here we we wish to present the qualitative arguments in favor of a possibility of localized 
solutions. The plausible scenario for the discret spectrum of the pomeron is as follows: The ker- 
nel K, is proportional to the probability of radiation of perturbative gluons, see Eq. (^). When 
the quarks of the colour dipole are a distance r 3> R c apart, the nonperturbative colour fields 
stretched between the quarks may strongly modify the vacuum and suppress the perturbative 
gluonic fluctuations on the nonperturbative background in the vicinity of quarks. Because the 
nonperturbative background fields are strongly anisotropic, this suppression of the perturba- 
tive gluonic fluctuations also will be anisotripic. However, to a very crude approximation this 
suppression can be modelled by the decrease with increasing r of the effective (nonlocal) pertur- 
bative coupling as{r,p) and/or the increase of hg with increasing r, resulting in the decrease 
of the kernel K, with increasing r. In terms of the Schrodinger equation fl42|) this amounts to a 
rise of the potential U(r) towards large r. (One can draw a useful analogy with the asymptotic 



freedom decrease of cts(r) and the related rise of the V(z) Eq. ([44D in the Schrodinger equation 
(P^).) As a natural scale at which such a suppression of the perturbative fluctuations can take 
place one can think of the confinement radius R con f, which is of the order of the diameter of 
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hadrons, R CO nf ~ 1 — 2f. Evidently, the energy e of the lowest state will be higher than the bot- 
tom of the potential well, so that the intercept Aip of the rightmost singularity will be lowered 
compared to the case of the minimal regularization. If with the increasing r the potential U(r) 
flattens at the still negative value U(oo) = — A c , then the Schrodinger equation (|29|) will have 
the continuous spectrum starting with e = — A c and in the complex j-plane there will be a cut 
from j = 1 + A c to j = — oo and certain number of isolated poles to the right of the cut, at 
1 + A c < j < 1 + Aip. If [/(oo) > 0, then the cut in the complex j-plane will start at j = 1. 

Although the above considerations are very qualitative ones, we regard them as giving a sound 
correlation between the discret spectrum of the pomeron and the nonperturative suppression of 
the perturbative gluonic fluctuations in large colour dipoles. (We do not suggest a dynamical 
mechanism for such a suppression, though.) In section 5.3 we explained why we disagree with 
Lipatov's conclusion [3] that the running coupling does by itself generate the discret spectrum 
of the pomeron (see section 5.3). Ross and Hancock [11] try to eliminate the contribution of 
the infrared region by the fc 2 -cutoff in the asymmetric form of the BFKL equation (p2|) and find 
a discrete spectrum of the pomeron for few models of the nonperturbative gluon propagator. 
The significance of their findings is not clear to us, since the Hancock-Ross procedure violates 
the symmetry of the cutoff of gluon propagators and its consistency with gauge invariance is 
questionable (for the criticism of this procedure see section 3). For this reason the correspondance 
between the Ross-Hancock cutoff and the above nonperturbative suppression of the perturbative 
gluon fluctuations is not evident to us. 

6 The evaluation of Aip and the pomeron dipole cross 
section 

The analytic solution of Eqs. ( |10| , |il"l) at finite R c and with the running coupling is not available, 
and we resort to the numerical analysis. We study the large-£ behavior of numerical solutions of 
Eq. (|IID and verify that cr(£, r) has the same asymptotic behavior cr(£, r) ==>- crjp(r) exp(Aip£) 
(within the overall normalization) irrespective of the boundary condition at £ = 0. Namely, we 
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compute the effective intercept 

A e// (£,r) = dlog<7(£,r)/0£ (52) 

and check that at large £ the effective intercept A e ff(£, r) tends to the same limiting value Ajp 
for all r. Flattening of A e jy(£,r) as a function of r and as a function of £ takes place at about 
the same value of £, and in all the cases the flattening was good to a few per mille. In Fig.l we 
present the eigenfuction ajp(r) for few values of the gluon correlation radius R c — 1/ Hg (we keep 
the same running strong coupling as(r) which freezes at oig = 0.8). As we have anticipated 
above, an»(r) flattens at large r and decreases towards small r faster than the BFKL solution 
flUf ), but slower that oc r 2 . 

The physical interaction picture behind this emergence of the factorizing r and £ dependences 
of the asymptotic cross section is as follows: Consider, for instance, interaction of two small-sized 
dipoles ta,tb <C R c - Radiation of perturbative gluons by such dipoles is suppressed by colour 
cancellations and by the small coupling ats(r). However, the perturbative gluons start sticking 
out of the initially small dipoles at a distance r ~ R c . Consequently, after few steps of radiation 
and/or after few units in £ the gluonic cloud builds up which has the size ~ R c , the evolution of 
which is controlled by as(R c ) ~ a si an d whose interaction properties do not depend on the 
size of the initial beam and target dipoles. The size of the beam and target dipoles only controls 
the rate of the evolution of the limiting gluonic cloud (for the more discussion on this point see 
below). 

In Fig. 2 we present the intercept Aip for few values of \iq = 0.3, 0.5, 0.75, 1.0 GeV. As we 
have emphasized above, this intercept is controlled by interactions of the semiperturbative gluons 
at transverse distances r ~ R c , where the strong coupling is rather large, a s (R c ) - af\ For 
the comparison, with as = 0.8 the BFKL formula ([15]) gives Ajp(BFKL) = 2.12. We remind 
that the BFKL formula is only valid in the scaling limit of R c — > oo and obviously overestimates 
the intercept for the case of the finite R c . In the literature one often cites the CK lower bound 
[9] 

Ajp > — a { s fr) (53) 

7T 

This bound was derived in the model with the CK sharp cutoff k 2 > in the asymmetric form 
of the BFKL equation (|22| ) with massless gluons, m 2 = 0. For a closer comparison with the CK 
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bound we substitute ctg in Eq. (|53| ) for the min {cts(fio), 0.8}. The resulting bound is shown 
in Fig. 2 and is not supported by our results. The significance of the CK bound is questionable, 
though, since the derivation [9] of the CK bound is plagued by violation of the symmetry of the 
BFKL equation, for a criticism see section 3. 

In terms of the exponent 7 of the small-r solutions (f48|), we have 7 — 1 = 1.19, 1.81, 2.33, 2.70 
at fie = 0.3, 0.5, 0.75, 1.0 GeV, respectively. The smaller is fiQ, the smaller is 7 — 1 and the 
steeper is the decrease of o~ip(r) at small r, in agreement with the results shown in Fig.l (for 
more discussion on the relevance of solutions (^) to o"ip(r) see below, section 7.3). 

In section 5.5 we discussed a possibility of generation of the discret spectrum of the pomeron 
via the (purely hypotetical) mechanism of the nonperturbative suppression of the perturbative 
gluonic fluctuations in colour dipoles of large size. Such a suppresion must lower the intercept 
of the pomeron. For a crude evaluation of the possible effect, we introduce in front of the kernel 
/C the suppression factor 



We consider the case of /ig = 0.75GeV, i.e., R c = 0.275f, and take d = 0.5f and R con f = 1.5f, 
i.e., impose strong suppression of the perturbative gluonic fluctuations already at r of the order 



factor corresponds to an abrupt rise of the potential U(r) at r « R con f. All this serves to enhance 
the discussed effect. We find a negligible, less than one per cent change of Aip with respect to 
the case of R CO nf = 00. Such a small effect could have easily been anticipated: the intercept Ajp 
is predominantly sensitive to the behavior of cr-jp(r) at r ~ R c , and although we have enforced 
very dramatic suppression of the perturbative gluonic fluctuations, the effect on the intercept is 
small because of the strong inequality R c <C R CO nf- This interpretation is confirmed by Fig. 3, 
in which we show the ratio of the (Jjp(R con f = 1.5/, r) to the crjp(R C onf = 00, r), and this ratio 
is essentially identical to unity up to r ~ 5R C ~ 1.4f. This result strongly suggests, that the 
intercept of the pomeron is the semiperturbative quantity, insensitive to the truly large-distance 
effects. 

In this paper we only concentrate on the rightmost singularity, and we do not evaluate the gap 
between the poles in the discret spectrum. We only comment, that the convergence of A e ff(£, r) 




2 



(54) 



1 + exp [0( r _ R conf )( r - R conf ) 2 /d 2 ] 



of the diameter of hadrons. In terms of the Schrodinger equation ( [4*2]) such a steep suppression 
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to the limiting value Aip becomes faster when the large-r cutoff (|54| ) is imposed, which suggests 
that the rightmost singularity is a pole separated from other singularities by a gap. 

Consequently, with all the due reservations on the large-distance QCD, we conclude that the 
intercept of the perturbative QCD pomeron is significantly higher than the effective intercept 
A ip (hN) ~ 0.1 as given by the phenomenology of hadronic scattering [15]. The plausible 
scenario, suggested by the observed slow rise of the hadronic cross sections, is that the exchange 
by perturbative gluons is only a small part of <Jtot{hN) at moderate energies ([5], for the early 
discussion of this scenario see [16]). The detailed phenomenology of the hadronic cross sections 
will be presented elsewhere. 



7 The GLDAP and BFKL evolutions in deep inelastic 
scattering 

7.1 Deep inelastic scattering and the dipole cross section 

The representation (3) is universal and also applies to the deep inelastic scattering [4,6]. The 
wave functions of the (T) transverse and (L) longitudinal virtual photon of virtuality Q 2 were 
derived in [6] and read 

|VMz ' r)|2 = J^f 2 ? Z2fi[z2 + zf]£2Kl{£r)2 + m2 f K ^ £r ^ > ( 55 ) 

R N f 

\*dz,r)\* = jglE^/ 2 Q 2 z\l-zfK,{erf , (56) 

where 

e 2 = «(1 - z)Q 2 + m) . (57) 

In Eqs. (|55|) -(|57D m./ is the quark mass and z is the Sudakov variable, i.e. the fraction of photon's 
light-cone momentum g_ carried by one of the quarks of the pair (0 < z < 1). Then 

a Tii ( 7 *iV,£,Q 2 ) = £dz J d 2 r\y TtL (z,r)\ 2 a(Z,r) (58) 

and the structure function is calculated as 

F 2 (^Q 2 ) = -^\a T ( 7 *N^,Q 2 ) + ( x i ( 7 *iV,e,Q 2 )l . (59) 



An 2 a e 
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Making use of the properties of the modified Bessel functions, after the ^-integration one can 
write 

M7*iV,£,Q 2 ) =! X dzf d*r\* T (z,r)\ 2 a(t,r) oc -L ^l^A ( 60 ) 

Jo J Q z Ji/Q 2 r z r l 

Notice, that the factor 1/Q 2 in Eq. (0), which provides the Bjorken scaling, comes from the 
probability of having the qq fluctuation of the highly virtual photon. The ratio cr(£,r)/r 2 is a 
relatively smooth function of r, which slowly rises towards small r, for the conventional DLLA 
behavior of cr(£, r) see Eqs. (|25| , |28|) , for the weak-coupling solution of our generalized BFKL 
equation see Eq. (|^). How the dipole cross section <r(£, r) is probed in deep inelastic scattering, 
will be discussed elsewhere ([17], see also [18]). Therefore, for the semiquantitative understanding 
of the transition from the GLDAP to the BFKL evolution we can concentrate on the r and £ 
dependence of cr(£,r). 



7.2 The transition between the BFKL and DDLA regimes 



One usually discusses the low-x behavior of structure functions in terms of the DLLA solution 
of the GLDAP evolution equations. In section 4 we already commented on the breaking of 
DLLA at large £ = log(^). Here we wish to concentrate on a comparison of the ^-dependence 
of the perturbative GLDAP and BFKL solutions for cr(£,r). The question of whether there is 
a strong, experimentally observable, difference between the BFKL and GLDAP evolutions, is 
being discussed in the literature for quite a time ([19] and references therein). 
The DLLA solution to GLDAP evolution (128]) gives the effective intercept 



A 



DDLA 



d_ 



4 



L(R,r) 



1 1 



(61) 



which must hold at 7] > 1 and moderately large £ > 1. The DLLA intercept (|6T|) vanishes at large 



£. On the other hand, in the BFKL regime the effective intercept fl52|) tends to Ajp irrespective 
of the radius r. This suggests the practical criterion: the DLLA approximation breaks at such 
values of £ = £ c ( r ), when A DLL a(£,,t) becomes smaller than Aip: 



A DLLA (£ c (r),r) < Ajp 



(62) 
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Notice, that since in expansion fl28|) the saddle point corresponds ton~ y/rj, this criterion when 



combined with the effective intercept essentially coincides with the DLLA breaking estimate 

For a meaningfull comparison of the GLDAP and BFKL evolutions, we must take the identical 
initial conditions. For the sake of definitness we take R c = 0.275f and consider scattering on the 
proton target, so that our initial condition for Eq. (|TTJ) is given by 

<t = m = |y rf 2 ^!*^^)! 2 ^^^) , (63) 

where Ru is a separation of quarks "1" and "2" in the proton (for the details see [6]). Since the 
DLLA asymptotics (^) works at £ > 1, we shall consider the cross section fl63|) as a result of 
the GLDAP evolution by £0 units from a lower energy, i.e., we take for the DLLA solution 



(Tdlla{^t) = a(0,r)J exp 



2J^h(R f ,r) (y'Cn -i ( ^ 



(64) 



It satisfies <Jdlla(£, — 0, r) = cx(0,r) by the construction, and we use it for evaluation of the 
DLLA intercept in Eq, fl5"T|). The same cr (r) is taken for the boundary condition for the BFKL 
equation (|TTJ). We assume this boundary condition to correspond to x = xq ~ 3 • 10~ 2 (the more 
detailed phenomenology of the deep inelastic scattering will be presented elsewhere [17]). 

Notice, that Eq. ( p5|) is the leading-logarithm formula defines L(Rf,r) up to an additive 
constant c < 1, which depends on the detailed form of the initial condition for ther GLDAP 
evolution, but which can be neglected in the DLLA limit of L(Rf, r) ^> 1. We compute A e yj(£ = 
0, r) from our BFKL equation (|Tl~D and make the readjustment 

ocs(r) 

such that Adlla{£, = 0, r) of Eq. (|6l|) gives a good approximation of A e //(^ = 0,r) at small r. 
With £0 = 1-25 this is achieved taking c ~ 0.05. 

The results of such a comparison of the BFKL and GLDAP evolutions are shown in Fig. 4. 
At £ = log(^r) ~ 1, the both DLLA and BFKL effective intercepts are smaller than Aip at 
r > 0.2f and are larger than Aip at smaller r. A good agreement of the BFKL and DLLA 
effective intercepts at small r is not surprizing, since our generalized BFKL equation (|11|) has 
the GLDAP equation as a limiting case at small r (also see a discussion in [5]). The ^-evolution 
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L(R f ,r) — > log 



+ c (65) 



of the DLL A and BFKL effective intercepts is very much different, though: The BFKL effective 
intercept starts flattening as a function of r and tends to Aip, rising at large r and decreasing 
at small r. In the opposite to that, the DLLA intercept monotonically decreases with £ at all r, 
until the GLDAP breaking fl62|) takes place. 

The boundary between the BFKL and DLLA regimes is given by our Eq. fl32|) . The intercept 
Ajp is numerically small, Ajp = 0.4 at \lq = 0.75, and the large numerical factor 

ik" s (66) 

emerges in the r.h.s. of Eq. (|32|) . Consequently, the DLLA to GLDAP evolution, and the 
GLDAP evolution by itself, may remain numerically viable in quite a broad range of £, relevant 
to the kinematical range of HERA. 

Closer inspection of Fig. 4 reveals a certain regularity: AuLLA{C,r) decreases with £ faster 
than A e ff(^,r), and the rate of the divergence is higher the larger is the radius r. Hence, we 
look at the combined r and £ dependence of cr(£,r). 

7.3 Testing the DLLA identity for the gluon distribution 

In the DLLA, the gluon structure function G(£, r) satisfies the equation (we consider Nf = 3 
active flavours) 

re \ 3 1 d 2 G(Z,r) 
K & r) = TG&r)dtdL(R,r) =1 - (6?) 

We shall refer to the equality /c(£, r) = I as the DLLA identity. One can easily evaluate «;(£, r) 
for the experimentally measured gluon distributions, and it is tempting to consider the departure 
from the DLLA identity as a measure of the accuracy of the GLDAP evolution. In this section 
we apply such a test of the DLLA identity to the above described solution of our generalized 
BFKL equation (|TT|) subject to the boundary condition (|63|) . 

The results of such a test are shown in Fig.5. We find that our BFKL solution produces 
r) ~ I in a very broad range of £ and r of the practical interest. Remarkably, the DLLA 
identity holds to the (20 — 30)% accuracy even at large r, up to r 2 < \R 2 C1 and to the few per 
cent accuracy at r < \R C - Notice somewhat oscillatory r-dependence of the A e //(£, r) for our 
BFKL solution. These oscillations die out at large £. They have an origin in the presence of 
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harmonics with large \v\ and/or large \/3\ in the expansion of the boundary condition ( |S^ ) in 
terms of eigenfunctions of the kernel /C (for instance, see Eqs. (|16| . |T7| ). These eigenfunctions give 
the oscillating contribution to cr(£, r), see Eqs. ( |ll3| , [40"|) . We have checked that the DLLA identity 



k(£, r) = 1 is satisfied by solutions of our generalized BFKL equation to a very good accuracy 
in a broad range of R c and Aq C d- 



The obvious conclusion from this observation is that the DLLA solution (2S), which is only 
valid at moderate values of £, evolves at larger £ into our new solution fl5Cf ). These new solutions 
satisfy our BFKL equation to the accuracy ~ as(r) and also satisfy the DLLA identity, although 
it was not obvious that necessarily ajp(£,,r) belongs to this new family of GLDAP/BFKL solu- 
tions (pUp. The exponent 7 is uniquely fixed by the pomeron intercept Ajp via Eq. (f£9p. We 
have checked that the pomeron dipole cross sections ajp(r) shown in Fig. 1 indeed to a good 
accuracy satisfy the property 

^M« s ( r )7-i = const ( 68 ) 

up to r < \R C . 

This observation shows that, as matter of fact, there is no real clash between the GLDAP and 
BFKL evolutions if as{r) is running. To his end, an important notice is that the discussion of the 
BFKL effects in the current literature concentrates upon the approximation of fixed as (for the 
review and references see [20]). We have found that the effect of the running coupling constant 



is quite substantial. Firstly, the pomeron cross section (50) dramatically differs from the oc r 1 
BFKL scaling solution at fixed as, see also Fig. 1. Secondly, one can define the counterpart of 
the DLLA identity for fixed as too, making the substitution 

Po f R2 dp 2 p (R 2 ' 

Then, the BFKL scaling pomeron solution fll8|) gives 



L(R, r) = £ / „ ^as(p) — ?Las log ^ . (69) 



K = 21og2, (70) 

which also emphasizes a dramatic difference between the cases of the fixed and running strong 
coupling. Therefore, matching the fixed-a^ BFKL regime at small x with the GLDAP regime 
at larger x is doomed to inconsistencies, and such considerations are not appropriate for the 
phenomenology of deep-inelastic scattering. 
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One interesting feature of the DLLA identity is worth while of mention. We have calculated 
ac(£, r) for solutions of our generalized BFKL equation for various boundary conditions, including 
the completely unphysical cr(£ = 0,r) which has the form of two Gaussians peaks with large 
spacing in r. Because of the diffusion in the radius r, the dip between the two Gaussians fills in 
very rapidly. As soon as such a <r(£, r) becomes a relatively smooth function of r, this is followed 
by a rapid approach of k(£, r) to unity. We find <10-15% departure of k(£, r) from unity already 
at £ >2-3 and < 5% departure at £ > 5. This fullfillment of the DLLA identity takes place 
when a(£,r) has the r-dependence still very different from the typical DLLA solution and/or 
the asymptotic pomeron cross section <Jjp(r). One plausible interpretation of this observation is 
that apart from the solution for the rightmost singularity at j = 1+Ajp, the (approximate) DLLA 
identity holds for other solutions in a relatively broad range of j < 1 + Ajp. The corresponding 
analysis goes beyond the scope of the present paper. 

We conclude this section by the prediction of the universal scaling violation at asymptotically 
large £. Namely, making use of (|50D in Eq. (|60D, we obtain 

F 2p (x,Q 2 ) oc 

(The possible difference between Q 2 of the virtual photon and the scale q 2 in the strong coupling 
in the r.h.s of this formula, and the conditions for the onset of this universal, ^-independent, 
scaling violation will be discussed elsewhere [17].) The results of Fig. 4 suggest, though, that 
the onset of this universal scaling violation is well beyond the kinematical range of the HERA 
experiments. 

8 The beam-target symmetry and the boundary condi- 
tion 

In the above discussion we have suppressed the target size variables r B but, evidently, it 
is the generalized energy- dependent dipole-dipole cross section <r(£, r^, r B ) which emerges as 
the fundamental quantitiy of the lightcone s-channel approach to the diffractive scattering. 
The lowest-order dipole-dipole cross section has an obvious beam-target symmetry property 



as(q 2 



'■An-. /]\ A IP 



(71) 
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o"o(^a,^*b) = &o{tb-iTa). This beam-target symmetry is but a requirement of the Lorentz- 
invariance and one has to have it extended to all energies. Thus, the beam-terget symmetry 
emerges as an important constraint on the admissible boundary condition for the BFKL equa- 
tion. In the derivation of our generalized BFKL equation we have treated the s-channel gluon 
g s of Fig. 6 as belonging to the beam-dipole A. The gluon-induced correction to the total 
cross section was reinterpreted in terms of the generalized dipole cross section cr(£, 7% r^), see 
Eq. (P). Alternatively, we could have treated the same s-channel gluon g s as belonging to the 
target-dipole, and the result must have been the same. The two descriptions differ in that in the 
former the perturbative t-cahnnel gluons g\ and gy enter the kernel /Q4, which we supply with 
the subscript A as it acts on the beam variable of the dipole-dipole cross section cr(£, r^, r B ), 
whereas the gluons gi and gy are the exchanged gluons in the dipole-dipole cross section ([[]). In 
the latter description #2, gy enter the kernel /Cg which now acts on the target variable tb of the 
dipole-dipole cross section cx(£, 7% r B ), and gi,gr will become the exchanged gluons in the input 
dipole-dipole cross section. The beam-target symmetry constraint essentially implies, that the 
boundary condition for the BFKL evolution must be calculable in the same perturbation theory 
as the one used to construct the generalized BFKL kernel K, Eq. (|T0|). (To this end we remind, 
that the Ki(x) in the kernel /C is precisely the derivative of the gluon propagator (correlation 
function) K (x).) This kernel-cross section relationship is crucial for having the beam-target 
symmetry, which will be violated with the arbitrary choice of the boundary condition for the 
dipole-dipole cross section cr(£ = 0,7% re). 



9 Restoration of the factorization at asymptotic energies 

It is well known that in the conventional Regge theory the isolated poles in the complex- j plane 
give rise to the factorizing scattering amplitudes. Even if the factorization holds for separate 
singularities in the j-plane, the scattering amplitude which is a sum of contributions from many 
close-by singularities, does not factorize. For instance, the lowest order dipole-dipole cross section 
([]]) manifestly breaks the factorization relation. It is interesting to notice, that the factorization 
of the total cross sections restores at asymptotic energies. 

Indeed, for the rightmost singularity in the j-plane the beam-target symmetric dipole-dipole 
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total cross section will have the factorized form 



cjp(£, r A ,r B ) = o-jp(r A )ajp(r B )exp(^Ajp) 



(72) 



By virtue of Eq. (0) this implies that at asymptotic energies 



a tot (AB) oc 



J dz A d 2 f A \^f(z A ,f A )\ 2 ajp(r A ) ■ J dz B d 2 f B \^/(z B ,f B )\ 2 ajp(r B 



exp(A IP 0, (73) 



which has the manifestly factorized form. From the point of view of the phenomenology of the 
total cross sections, this factorization will be broken by the unitarization corrections needed 
to tame too rapid a rise of the perturbative bare pomeron cross section. The above predicted 
universal scaling violation in deep inelastic scattering is a particular case of this factorization 
restoration. 



10 Measuring the BFKL pomeron at HERA 

The experimental determination of the pomeron intercept Ajp requires isolation of the BFKL 
pomeron cross section, which is not an easy task because of the close similarity of the BFKL and 
GLDAP solutions, discusssed in detail in Section 7. The GLDAP evolution is a special limiting 
case of the generalized BFKL evolution. Every and each solution of the BFKL equation does not 
satisfy the GLDAP evolution and the DLLA identity. Our remarkable finding is that the small- 



r pomeron cross section (|50|) which is a solution of the BFKL equation, satisfies the GLDAP 
evolution equation within the correction terms ~ cts(r). Consequently, one may expect that 
GLDAP evolution with the properly posed boundary condition can be a good approximation 
to the BFKL evolution. The GLDAP considerations can not by themselves fix this boundary 
condition, though. Inspection of Fig. 4 shows that the DLLA and BFKL effective intercepts are 
close to each other and to the Ajp in a broad range of x of the practical interest at 

r^r = (0.4 - 0.5) R c . (74) 

Consequently, choosing boundary condition G(£,r ) oc x~ Ajp at r »s r , one will obtain the 
GLDAP solutions, which at r < r$ for all the practical purposes will be indistinguishable from 
the BFKL solution. Eq.(|K]) shows that the proton structure function receives contributions from 
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the broad range of r, in which the effective intercept A e //(r) is substantially different from Aip. 
For this reason studying the x-dependence of Fz(x, Q 2 ) is not suitable for accurate determination 
of the pomeron intercept. 

Fig.4 suggests that zooming at r ps0.15-0.2f one can determine the pomeron intercept Ajp 
already at the moderately large - of the HERA experiments. This can be achieved measuring 
either the longitudinal structure function Fl(x,Q 2 ) or the scaling violations in the transverse 
structure function dF T (x, Q 2 )/d\og(Q 2 ). Making use of the wave functions (|55| , |56|) one can easily 
show that these quantities are local probes of <r(£, r) at 

r 2 ~ ^ , (75) 

where for the longitudinal structure function Bl ~ 12 and Bt ~ 5 for the transverse structure 
function. Consequently, experimentally Ajp can best be estimated as 

d\ogF L (x,Q 2 ) 

Ajp = — (76) 

a log a; 

at Q 2 = Ql ^12-25 GeV 2 and/or 

A 1 9 2 log F T (x,Q 2 ) 

w F T (x,Q 2 ) d log x d log Q 2 1 ' 

at Q 2 = Q\ ~5-8GeV 2 . Our prediction is that these derivatives must be x- independent already 
starting with x < 10~ 2 . The above difference by the factor w (2.5 — 3) between Q\ and Q\ is an 
important consequence of wave functions (|55|j56|) . Purely kinematically, the smaller Q 2 implies 
the broader range of -, and for this reason the second determination (|77|) may prove more 
accurate one. A comparison of the two determinations ([76]) and ( [771) is an important consistency 
check. We conclude that from the point of view of measuring the pomeron intercept, the HERA 
experiments must concentrate on accurate measurements of F 2 (x, Q 2 ) and Fl(x, Q 2 ) in the region 
of moderate Q 2 < 30 GeV 2 . 



11 Conclusions 



The purpose of this paper has been to understand the spectrum of eigenvalues of the generalized 
BFKL equation [4,5] for the dipole total cross section. Our emphasis was on the realistic case of 
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the finite correlation radius for the perturbative gluons and of the freezing strong coupling. The 
advantage of our generalized BFKL equation ( |T0| , [TT| ) is an easy introduction of the finite gluon 
correlation radius in the manner which is consistent with the gauge invariance constraints. 



We have shown that our generalized BFKL equation (|10|,|TT]) has the continuous spectrum, 
which corresponds to the QCD pomeron described by the cut in the complex j-plane. We have 
determined a dependence on the size r of the dipole for the rightmost singularity in the j-plane 
and found the corresponding intercept An>. It is much smaller than given by the BFKL formula 
fllED and even smaller than the lower bound cited in [9], but is substantially larger than the 
phenomenological value Ajp(hN) ~ 0.1. we have investigated a transition from the GLDAP to 
the BFKL evolution in the limit of large -, and have shown that the two equations have an 
overlapping asymptotic solution in the weak coupling limit. The phenomenologically important 
implication is that the GLDAP evolution with the correctly posed boundary condition remains 
a viable description of deep inelastic scattering well beyond the kinematical range of the HERA 
experiments. Nonetheless, the pomeron intercept can be determined already from the HERA 



data on deep inlastic scattering, and we have suggested practical methods of measuring A 



ip- 
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Figure captions: 



Fig.l - The pomeron dipole cross section cxip(r) for different values of \xq- The straight lines show 
the r 1 and r 2 behavior. 

Fig. 2 - The intercept An> for \iq = 0.3, 0.5, 0.75, 1.0 GeV (shown by triangles). The solid curve 
shows the Collins-Kwiecinski lower bound (|53|). 



Fig. 3 - The ratio of the pomeron dipole cross section with the cutoff (|54T ) of the large-r controbu- 
tion to the pomeron dipole cross section without the cutoff. The cutoff radius R con f = l-5f, 
the gluon correlation radius R c = 0.275f. 



Fig. 4 - Comparison of effective intercepts of the DLLA solution (|64l ) to the GLDAP evolution and 
of the solution of our BFKL equation. Both solutions start with the identical dipole cross 
section at x = 3 • 10 -2 . The pomeron intercept Ajp = 0.4 is shown by the horizontal line. 

Fig. 5 - Test of the DLLA identity for the solution of our generalized BFKL equation. The gluon 
correlation radius R c = 0.2751 

Fig. 6 - One of diagrams of the driving term of the rising dipole cross section. 
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